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ABSTRACT 

In the Local Group, nearly all of the dwarf galaxies (-Mstar ^ 10® M©) that are satellites within 
300 kpc (the virial radius) of the Milky Way (MW) and Andromeda (MSI) have quiescent star for¬ 
mation and little-to-no cold gas. This contrasts strongly with comparatively isolated dwarf galaxies, 
which are almost all actively star-forming and gas-rich. This near dichotomy implies a rapid trans¬ 
formation of satellite dwarf galaxies after falling into the halos of the MW or MSI. We combine 
the observed quiescent fractions for satellites of the MW and MSI with the infall times of satellites 
from the Exploring the Local Volume in Simulations (ELVIS) suite of cosmological zoom-in simula¬ 
tions to determine the typical timescales over which environmental processes within the MW/M31 
halos remove gas and quench star formation in low-mass satellite galaxies. The quenching timescales 
for satellites with Mgtar < 10® Mq are short, < 2Gyr, and quenching is more rapid at lower Mgtar- 
These satellite quenching timescales can be 1 — 2 Gyr longer if one includes the time that satellites 
were environmentally preprocessed by low-mass groups prior to MW/M31 infall. We compare with 
quenching timescales for more massive satellites from previous works to synthesize the nature of satel¬ 
lite galaxy quenching across the observable range of Mgtar = 10®“^^ Mq. The satellite quenching 
timescale increases rapidly with satellite Mgtar, peaking at « 9.5 Gyr for Mgtar ~ 10® Mq, and the 
timescale rapidly decreases at higher Mgtar to < 5 Gyr at Mstar > 5 x 10® Mq. Overall, galaxies with 
Mstar ~ 10® Mq, similar to the Magellanic Glouds, exhibit the longest quenching timescales, regardless 
of environmental or internal mechanisms. 

Subject headings: galaxies: dwarf — galaxies: evolution — galaxies: groups: general — galaxies: star 
formation — Local Group — methods: numerical 


1. INTRODUCTION 

Galaxies in denser environments are more likely to 
have suppressed (quiescent) star formation and little-to- 
no cold gas than galaxies of similar stellar mass, Mgtar, 
in less dense environments. The observed environmen¬ 
tal effects within the Local Group (LG), on the satellite 
galaxies within the halos of the Milky Way (MW) and 
Andromeda (M31), are particularly strong (e.g., Einasto 
et al. 1974; Grcevich & Putman 2009; McConnachie 2012; 
Phillips et al. 2014; Slater & Bell 2014), even compared 
to the already strong effects on (more massive) satel¬ 
lites within massive groups/clusters (e.g., Wetzel et al. 
2012). Specifically, dwarf galaxies around the MW/M31 
show a strikingly sharp and nearly complete transition 
in their properties within « 300 kpc (approximately the 
virial radius, Rvu, of the MW or MSI), from irregular 
to spheroidal morphologies, from significant to little-to- 
no cold atomic gas, and from star-forming to quiescent. 
This trend has just a few exceptions: 4 gas-rich, star¬ 
forming galaxies persist within the halos of the MW (the 
LMC and SMG) and MSI (LGS 3 and IG 10), and 4-5 
quiescent, gas-poor galaxies reside well beyond i?vir of ei¬ 
ther the MW or MSI: Getus (Lewis et al. 2007), Tucana 
(Fraternali et al. 2009), KKR 25 (Makarov et al. 2012), 
KKs 3 (Karachentsev et al. 2015), and possibly An- 
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dromeda XVHl, though Getus and Tucana may have or¬ 
bited within the MW halo (Teyssier et al. 2012/ This ef¬ 
ficient satellite quenching is particularly striking because, 
other than KKR 25 and KKs 3, at Mgtar < 10® Mq all 
known galaxies that are sufficiently isolated (> 1500 kpc 
from a more massive galaxy) are star-forming (Geha 
et al. 2012; Phillips et al. 2014). Thus, the MW and 
MSI halos show the strongest environmental influence 
over their satellites of any known systems, making the 
LG a compelling laboratory for studying environmental 
processes on galaxies. 

Several such processes within a host halo can regulate 
the gas content, star formation, morphology, and even¬ 
tual disruption of satellites, including gravitational tidal 
forces (e.g., Dekel et al. 2003), galaxy-galaxy tidal inter¬ 
actions (e.g., Farouki & Shapiro 1981), galaxy-galaxy 
mergers (e.g., Deason et al. 2014), and ram-pressure 
stripping of extended gas (e.g., McGarthy et al. 2008) 
or inter-stellar medium (e.g., Gunn & Gott 1972; Ton- 
nesen & Bryan 2009). The key astrophysical challenge is 
understanding the relative importance of these, includ¬ 
ing which (if any) dominate, and how they vary across 
both satellite and host masses. 

One strong constraint comes from determining the 
timescale over which environmental quenching occurs, as 
previous works explored at higher masses (e.g., Balogh 
et al. 2000; De Lucia et al. 2012; Wetzel et al. 2013; 
Hirschmann et al. 2014; Wheeler et al. 2014). For the 
satellite dwarf galaxies in the LG, recent works showed 
that their environmental quenching efficiency is higher 
than for higher-mass satellites (Phillips et al. 2014; Slater 
& Bell 2014). In this letter, we combine the observed 
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Figure 1. For all known dwarf galaxies in the Local Group out 
to 1.6 Mpc, the distance from their nearest host (MW or M31) 
versus stellar mass, Mgtar? or absolute magnitude, My. Points 
show individual galaxies: actively star-forming (blue stars) and 
quiescent (Afgas/Afstar < 0.1, red circles). Black curve shows the 
detection limit for dwarf spheroidal-like galaxies with stars resolved 
in SDSS, plus extrapolation beyond ~ 700 kpc. 

quiescent fractions for satellites in the MW/M31 halos 
with their typical infall times from cosmological simu¬ 
lations to infer the timescales over which environmental 
processes remove their gas and quench star formation. 
Motivated by Wetzel et al. (2015), we also consider the 
possible impact of group preprocessing on satellites be¬ 
fore they fell into the MW/M31 halos. We also compare 
with previous works on more massive satellites, to syn¬ 
thesize satellite quenching across the observable range of 

2. METHODS 
2 . 1 . Observations 

To examine the observed properties of dwarf galaxies 
in the LG, we use the compilation from McConnachie 
( 2012 ), which includes all galaxies known at that time 
within 3 Mpc of the Sun. We also include the more recent 
observations of cold atomic gas mass from Spekkens et al. 
(2014). We define “satellite” galaxies as those within 
300 kpc of either the MW or M31, motivated by the ob¬ 
served sharp transition in star formation, gas mass, and 
morphology within this distance. 

Observed dwarf galaxies show a tight correlation be¬ 
tween morphology, star formation, and cold gas mass: 
all spheroidals have little-to-no detectable cold gas (e.g., 
Spekkens et al. 2014) or star formation (e.g., Weisz 
et al. 2014a), and almost all irregulars have significant 
cold gas and ongoing star formation. Thus, we de¬ 
fine “quiescent” galaxies as having Mgas/Mgtar < 0.1 
or, if they have no cold gas constraints, having col¬ 
ors/morphologies that resemble spheroidals. By this 
definition, the only star-forming, gas-rich satellites are: 
LMC (Mstar = 1-5 X 10® M 0 , Mgas/Mstar « 0.3) and 
SMC (Mstar = 4.6 X 10® Mq, Mgas/Mstar ^ 1) around the 
MW, LGS 3 (Mstar = 9.6 X 10® M©, Mgas/M,tar « 0.4) 
and IC 10 (Mgtar = 9 x 10^, Mgas/Mgtar ~ 0.6) around 
M31. 

For each dwarf galaxy out to 1.6 Mpc, Figure 1 shows 


its distance from nearest host (MW or M31) versus Mgtar- 
Almost all quiescent galaxies are within « 300 kpc of 
their host. The black curve shows the detection limit 
(and extrapolation) for dwarf spheroidal-like galaxies in 
SDSS (Tollerud et al. 2008), which highlights complete¬ 
ness at different Mgtar- 

2 . 2 . Simulations 

To measure the infall times of satellites, we use Ex¬ 
ploring the Local Volume in Simulations (ELVIS), a suite 
of cosmological zoom-in V-body simulations intended to 
model the LG (Garrison-Kimmel et al. 2014) in ACDM 
cosmology: as = 0.801, ^matter = 0.266, Da = 0.734, 
rig = 0.963 and h = 0.71. Within the zoom-in regions, 
the particle mass is 1.9 x 10® M© and the Plummer- 
equivalent force softening is 140 pc physical. 

ELVIS contains 48 dark-matter halos of masses similar 
to the MW or M31 (Mvir = 1-0 — 2.8 x 10^^ M©), with 
a median i?vir ~ 300 kpc. Half of the halos are in a pair 
that resemble the masses, distance, and relative velocity 
of the MW-M31 pair, while the other half are single iso¬ 
lated halos. Given the lack of systematic differences in 
satellite infall times for the paired versus isolated halos 
(Wetzel et al. 2015), we use all 48 to improve statistics. 

ELVIS identifies dark-matter (sub)halos using the 
six-dimensional halo finder rockstar (Behroozi et al. 
2013b). For each halo, we assign a virial mass, Mvir, 
and radius, i?vir, according to Bryan & Norman (1998). 
A “subhalo” is a halo whose center is inside i?vir of a 
more massive host halo, and a subhalo experiences “first 
infall” and becomes a “satellite” when it first passes 
within i?vir- For each subhalo, we compute the peak 
mass, Mpeak, that it ever reached, and we assign M^tar 
to subhalos based on Mpeak using the relation from abun¬ 
dance matching in Garrison-Kimmel et al. (2014), which 
reproduces the observed mass function in the LG if one 
accounts for observational incompleteness (Tollerud et al. 
2008; Hargis et al. 2014). 

For more on ELVIS and its satellites’ infall times, see 
Garrison-Kimmel et al. (2014) and Wetzel et al. (2015). 

3. RESULTS 

3.1. Observed Quiescent Fractions for Satellites 

Eigure 2 shows, for all satellite galaxies at Mgtar ^ 
10® M© within 300 kpc of the MW or M31, the fraction 
that are quiescent, in 1-dex bins of Mgtar (see also Phillips 
et al. 2014; Slater & Bell 2014). We do not correct for 
observational completeness versus Mgtar (Figure 1), be¬ 
cause we measure the relative fraction in each bin, which 
is likely unbiased. We show fractions for all satellites 
(blue circles) and separately for those in the MW (vio¬ 
let squares) and M31 (green triangles) halos. Error bars 
show 68 % uncertainty for the binomial counts using a 
beta distribution. Of the 56 satellites, only 4 (7%) are 
star-forming/gas-rich: LMC and SMC of the MW, LGS 
3 and IC 10 of M31. Moreover, at Mgtar < 8 x 10^ M©, 
only 1 (LGS 3) of the 51 satellites is star-forming, and 
at Mgtar < 9 X 10® M© all 40 satellites are quiescent. 

These near-unity quiescent fractions for satellites of 
the MW/M31 contrast strongly with the nearly zero 
quiescent fraction for isolated (non-satellite) galaxies at 
Mgtar < 10® M© (Geha et al. 2012; Phillips et al. 2014). 
The only clear exceptions are the quiescent galaxies 
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Figure 2. For all satellites galaxies with Mstar $ 10® Mq within 
300 kpc of the Milky Way (MW) or Andromeda (M31), the fraction 
that are quiescent (Mgas/Afgtar < 0.1) versus stellar mass, Mgtar- 
Blue circles show all satellites, violet squares (green triangles) show 
those of just the MW (M31). Of these 56 satellites, only 4 are star¬ 
forming/gas-rich: LMC (Mstar = 1.5x 10® M0) and SMC (Mstar = 
4.6 X 10® Mq) around the MW, LGS 3 (Mstar = 9.6 x 10^ Mq) and 
IC 10 (Mstar = 9 X 10^ Mq) around M31. At Mstar < 8 x lO"^ Mq, 
50 of 51 satellites are quiescent, and at Mstar < 9 x 10^ Mq all are 
quiescent. Error bars show 68% uncertainty from observed counts. 

KKR 25 (Mstar = 1-4 x 10® M©) and KKs 3 (M^tar = 
2.3 X 10^ Mq) at « 2 Mpc from the MW/M31. (Though, 
as Figure 1 shows, the completeness distances at low 
Mstar leave open the possibility for more isolated qui¬ 
escent dwarf galaxies.) 

3.2. Inferred Quenching Timescales for Satellites 

We now translate the quiescent fractions in Figure 2 
into the typical timescales over which environmental pro¬ 
cesses quenched satellites of the MW/M31 after they fell 
into a host halo, following the methodology of Wetzel 
et al. (2013). 

First, motivated by the dearth of isolated galaxies with 
Mstar < 10® Mq that are quiescent at z « 0, we as¬ 
sume that all satellites with Mstar (z = 0) < 10® Mq were 
star-forming prior to first infall. However, we do not 
model Mstar (^ = 0) < IO'^Mq, because cosmic reion¬ 
ization likely quenched most/all such galaxies at high 
redshift (e.g., Weisz et al. 2014a; Brown et al. 2014). At 
Mstar(x = 0) = 10^“® Mq, Satellites’ star-formation his¬ 
tories show a mix of complete quenching by z > 3 (e.g., 
Bootes I, Leo IV) and signs of star formation at z < 1 
(e.g.. And XI, And XII, And XVI) (Weisz et al. 2014a,b; 
Brown et al. 2014), so quenching at these masses may 
arise from a mix of reionization and the host-halo en¬ 
vironment. That said, the 100% quiescent fraction for 
satellites at this M^tar means that if both processes are 
responsible, both are highly efficient. Furthermore, if the 
satellites that were quenched by reionization versus the 
host-halo environment have similar infall-time distribu¬ 
tions, our modeling approach remains valid. Thus, we 
include this Mstar but label it distinctly to emphasize 
caution in interpretation. 

Within each 1-dex bin of Mstar, we use ELVIS to com¬ 
pute the distribution of infall times for satellites at z = 0. 
Infall into the MW/M31 halo (or any host halo) typically 
occurred 5 —8 Gyr (or 7—10 Gyr) ago, and our most mas¬ 


sive satellites typically fell in 2—3 Gyr more recently than 
our least massive (see Figures 1 and 2 in Wetzel et al. 
2015). Assuming that environmental quenching corre¬ 
lates with time since infall, we designate those that fell 
in earliest as having quenched, and we adjust the time- 
since-infall threshold for quenching until we match the 
observed quiescent fraction in each Mgtar bin. 

Several works have shown that this model successfully 
describes the dependence of satellite quiescent fractions 
on host-centric distance (e.g., Wetzel et al. 2013, 2014; 
Wheeler et al. 2014) because infall time correlates with 
host-centric distance (e.g., Wetzel et al. 2015). However, 
this correlation means that we must account for the dis¬ 
tances of the observed satellites in computing their in¬ 
fall times. Thus, in ELVIS we only select satellites out 
to the maximum host-centric distance that they are ob¬ 
served in each Mgtar bin. This matters most at the high¬ 
est Mstar, where all observed satellites (M32, NGC 205, 
LMC/SMG) reside < 60 kpc from the MW or M31. 

Figure 3 shows the resultant environmental quenching 
timescales (the time duration from first infall to being 
fully quiescent/gas-poor) for satellites versus their Mstar 
(or subhalo Mpeak)- Blue circles show satellites in the 
MW and M31, and we shade the lowest Mstar to high¬ 
light caution in interpretation because of reionization. 
We derive error bars from the 68% uncertainty in the 
observed quiescent fractions in Figure 2. 

As explored in Wetzel et al. (2015), many satellites 
first fell into a another host halo (group), typically with 
Mvir 10^^ Mq, before falling into the MW/M31 ha¬ 
los. Such groups may correspond to, for example, the 
LMC, as the newly discovered dwarf galaxies near the 
LMC (Koposov et al. 2015; Bechtol et al. 2015) suggest, 
although such groups disperse in phase space ^ 5 Gyr 
after MW/M31 infall (Deason et al. 2015), and typically 
half of such preprocessing hosts do not survive to z = 0 
(Wetzel et al. 2015), so preprocessed satellites are not 
always easily distinguishable. Because the importance of 
this preprocessing in low-mass groups remains unclear, 
we present quenching timescales both neglecting (left 
panel) and including (right panel) group preprocessing. 
The latter results in longer quenching timescales, though 
it primarily shifts the upper tail of the distribution and 
not the median. 

Both panels show shorter median quenching timescales 
for less massive satellites: ^ 5 Gyr at Mgtar = 10®“® Mq, 
2 —3 Gyr at Mgtar = 10^“® Mq, and < 1.5 Gyr at Mgtar < 
10^ Mq, depending on group preprocessing. Moreover, 
the median timescale for two of the lowest Mgtar bins 
is 0 Gyr, because 100% of those satellites are quiescent, 
which implies extremely rapid quenching after infall. 

Figure 3 also shows infall/quenching timescales that 
are more directly measured for satellites of the MW. 
The 3-D orbital velocity measured for the LMC/SMG 
strongly suggests that they are experiencing first infall 
and crossed inside i?vir of the MW « 2 Gyr ago (Kallivay- 
alil et al. 2013). Given that the LMC and SMC remain 
star-forming, this places a lower limit to their quenching 
timescale (gray triangle), consistent with our statistical 
timescales. Similarly, measurements of the 3-D orbital 
velocity and star-formation history for Leo I indicate that 
it fell into the MW halo « 2.3 Gyr ago and quenched 
« 1 Gyr ago (near its « 90 kpc pericentric passage), im- 
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Figure 3. Satellite quenching timescales for galaxies across the observable range of stellar mass, Mgtar (top axis shows subhalo Mpeak from 
abundance matching). Blue circles show satellites of the MW and M31, obtained by matching the observed quiescent fractions in Figure 2 
to rank-ordered infall times of satellites from the ELVIS simulations (Wetzel et al. 2015) in 1-dex bins of Mgtar. At ALstar = Mq 

(light blue), reionization may have quenched some satellites prior to infall. Error bars come from the 68% uncertainty in observed quiescent 
fractions in Figure 2. Left panel uses time since first infall into the current MW/M31-like halo, while right panel uses time since first 
infall into any host halo, thereby including possible effects of group preprocessing. Gray triangle shows lower limit for the LMC/SMC 
system from its measured orbit (Kallivayalil et al. 2013), and gray pentagon shows the quenching timescale for Leo I from its measured 
orbit and star-formation history (Sohn et al. 2013). Red squares show times inferred for satellites with Mstar = 10®'®, 10®'® Mq around 
hosts with Matar > 2.5 X 10*® Mq in SDSS (Wheeler et al. 2014), and green curve shows the same for more massive satellites in groups 
of Mvir = 10*7“*® Mq in SDSS (Wetzel et al. 2013). The satellites in the MW/M31 halos quenched more rapidly after infall than more 
massive satellites (around other hosts). Overall, the quenching timescale increases with Mstar, is longest at Mstar ~ 10® Mq (near the 
masses of the Magellanic Clouds), then decreases with further increasing Mstar. 


plying a quenching timescale of « 1.3 Gyr (Sohn et al. 
2013, gray pentagon), again consistent with our results. 

The mass trend in Figure 3 is broadly consistent 
the star-formation-history-based results of Weisz et al. 
(2015) that more massive dwarf galaxies in the LG 
quenched more recently. Also, the overall timescale is 
broadly consistent with Slater & Bell (2014), who in¬ 
ferred a typical quenching time since first pericenter of 
1 — 2 Gyr, which implies a quenching time since infall of 
~ 3 Gyr, though they did not examine mass dependence. 

We also compare these timescales with previous re¬ 
sults for more massive satellites of other hosts. Fig¬ 
ure 3 (green curve) shows the quenching timescales for 
satellites in groups with Mvir = 10*2“*® Mq from Wetzel 
et al. (2013), who used identical methodology based on 
galaxies in SDSS (Tinker et al. 2011; Wetzel et al. 2012). 
(Hirschmann et al. (2014) found similar timescales ver¬ 
sus Mstar-) Red squares show timescales from Wheeler 
et al. (2014), who also used an SDSS galaxy catalog 
(Geha et al. 2012) and similar methodology for satellites 
with Mstar ~ io®'5“9-5 Mq around hosts with Mgtar > 
2.5 X 10*8 Mq, or Mvir ~ io*2'8“*“ Mq, much more mas¬ 
sive than the MW/M31. Both works measured satellite 
infall times, including group preprocessing, from cosmo¬ 
logical simulations. The timescale changes rapidly be¬ 
tween these works, from « 5.2 Gyr at Mgtar ~ lO®'® Mq 
to « 9.5 Gyr at « IO^ ^Mq. Both analyses used sim¬ 
ilar galaxy catalogs and methodologies, though Wetzel 
et al. (2013) used a group catalog to narrow the masses 
of the hosts, which are more similar to the MW/M31, 
while the hosts in Wheeler et al. (2014) are more mas¬ 
sive, on average. Thus, part of this change in timescale 


could arise if more massive hosts quench satellites with 
Mstar ^ 10® Mq less rapidly. Absent that, these results 
imply that the satellite quenching timescale rises rapidly 
near Mstar ~ 10® Mq. Furthermore, the timescale from 
Wheeler et al. (2014) implies some tension with our 
« 5 Gyr at Mstar ~ 10® ^Mq, as driven by the higher 
quiescent fraction in the MW/M31 at this Mstar, specifi¬ 
cally, NGC 205 and M32, two quiescent satellites of M31. 
This tension could be explained if NGC 205 and M32 
both fell into the M31 halo unusually early (> 9.5 Gyr 
ago), and/or (again) if M31 quenches its satellites more 
rapidly than the higher-mass (Mvir = 10*2'5-i4 Mq) host 
halos in Wheeler et al. (2014). 

Altogether, Figure 3 indicates a complex dependence of 
the satellite quenching timescale on Mstar- The typical 
timescale for satellites of the MW/M31 increases with 
-41star, from < 1 Gyr at Mstar < 10'’ Mq to ~ 5 Gyr 
at Mstar ~ 10®'5 Mq. Wheeler et al. (2014) indicates 
that this mass dependence continues, though with a 
rapid increase (~ 2x) to « 9.5 Gyr, and no change from 
Mstar « 10®-® to 109-5 Mq. Finally, Wetzel et al. (2013) 
shows that the timescale decreases near 5 x 10® Mq and 
continues to decline with increasing Mstar- Overall, the 
typical satellite quenching timescale is shortest at low¬ 
est Mstar, short at the highest Mstar, and longest at 
Mstar ~ 10® Mq, Comparable to the Magellanic Clouds. 

4. DISCUSSION 

We conclude by discussing the dependence of satellite 
quenching timescales on Mstar from Figure 3 in the con¬ 
text of the underlying physics. 

At Mstar ^ 10® Mq, the long timescales suggests that 
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satellite quenching is caused by gas depletion in the 
absence of cosmic accretion, via the stripping of ex¬ 
tended gas around a satellite after infall (“strangula¬ 
tion”). This scenario can explain shorter timescales at in¬ 
creasing Mstar, because higher-Mgtar star-forming galax¬ 
ies have lower Mgas/Mgtar (in cold atomic and molecu¬ 
lar gas, e.g., Huang et al. 2012; Boselli et al. 2014) and 
thus shorter gas depletion timescales in the absence of 
accretion. For example, Bradford et al. (2015) found 
that isolated galaxies follow Mgas/Mgtar cx Mg^ar^^ at 
Mstar > 10® ® M 0 . Furthermore, star-forming galaxies at 
ihfstar ^ 10® Mq have M^as = Mm + ~ i^star, with 

gas depletion timescales comparable to a Hubble time. 
Thus, satellite quenching timescales at Mgtar ^ 10® Mq 
do not necessarily require strong environmental processes 
beyond truncated gas accretion (see also Wetzel et al. 
2013; Wheeler et al. 2014; McGee et al. 2014). Further¬ 
more, at Mstar > 10 ® Mq, internal feedback from stars 
and/or black holes also may quench satellites after infall, 
which could help explain the shortening of the timescale 
with increasing M^tar- 

However, strangulation cannot explain the rollover in 
quenching times at Mgtar ^ 10® Mq, because the star¬ 
forming dwarf galaxies in the LG also have M^as ^ Mjtar 
(Grcevich & Putman 2009), enough to fuel star for¬ 
mation for a Hubble time. Thus, the rapid decline of 
the timescale at lower Mgtar requires an additional pro- 
cess(es) to remove gas from satellites after infall. This 
likely arises from increased efficiency of ram-pressure 
stripping in removing cold gas from such low-mass galax¬ 
ies, which have shallower potential wells. Moreover, the 
same internal stellar feedback that regulates the low star- 
formation efficiency in dwarf galaxies likely heats/drives 
significant cold gas to large radii (e.g., Muratov et al. 
2015), which would assist such environmental stripping. 
Thus, the rapid quenching timescales for dwarf galaxies 
may arise from the nonlinear interplay of both internal 
feedback and external stripping (e.g., Nichols & Bland- 
Hawthorn 2011; Bahe & McCarthy 2015). 

Overall, satellites with Mgtar ~ 10® Mq (similar to 
Magellanic Clouds) represent the transition between 
quenching via gas consumption and via gas stripping, 
and no quenching mechanism, either internal or exter¬ 
nal, appears to operate efficiently near this mass (see 
also Geha et al. 2012; Weisz et al. 2015). 

Finally, the above scenario may explain the curious 
similarity between the mass dependence of the quench¬ 
ing timescale in Figure 3 and the underlying galaxy-halo 
-^starZ-^haio ratio, which also is small at both high and 
low Mstar and peaks at Mstar 10^° Mq (e.g., Behroozi 
et al. 2013a). In particular, at high M^tar, the same 
physical process(es) that lowers Mstar/M^aio also lowers 
a galaxy’s cold gas mass, which in turn causes more mas¬ 
sive satellites to quench more rapidly, absent accretion. 
At low Mstar, the same shallower potential well that al¬ 
lows stellar feedback to lower Mstar/Mhaio also allows ex¬ 
ternal stripping to occur more easily and thus quenching 
to occur more rapidly. 

During preparation, we learned of Fillingham et al. 
(2015), who also used ELVIS to constrain the quenching 
timescales of satellites of the MW/M31 and reached 
similar conclusions. 
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